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INTRODUCTION 


The  use  of  an  organic  plasticizer  in  solid  rocket  propellant  binders, 
although  necessary,  has  occasionally  resulted  in  poor  rocket 
performance.  The  primary  reason  for  this  fault  was  the  loss  of 
plasticizer  from  within  the  fuel  pellet.  The  actual  mechanism  for 
this  loss  has  been  unknown.  Thus,  the  problem  of  plasticizer 
disappearance  within  a fuel  pellet  was  the  topic  of  this  research. 

Initially,  four  desired  results  were  envisioned  that  would  essentially 
complete  this  project.  These  objectives  were  as  follows; 

1.  Determination  of  factor(s)  resulting  in  plasticizer 
disappearance  (diffusion  or  kinetic  controlled  phenomena?) 

2.  Determination  of  physical  and/or  chemical  parameters 
descriptive  of  the  process  involved  in  the  migration/ 
reaction  of  said  plasticizer. 

3.  A mathematical  model  and  solution  for  the  concentration  of 
plasticizer  in  a fuel  pellet  as  a function  of  time  and 
whatever  special  dependence  was  appropriate. 

4.  Use  of  the  above  results  to  retard  or  prevent  the 
physical/chemical  process  from  occurring. 

While  some  of  these  results  would  undoubtedly  require  subjective 
evidence  to  support  various  conclusions  and  assumptions,  it  was  highly 
desirable  to  obtain  as  many  specific  formulations  and  results  as 
possible . 

DISCUSSION 

Fuel  Pellet  Construction/Characteristics 

.As  it  would  be  necessary  to  eventually  have  experimental  data  from  an 
actual  fuel  pellet  system  it  w'as  decided  to  fabricate  a pellet  of 
materials  similar  to  those  utilized  in  the  construction  of  the  fuel 
pellet  employed  in  the  Minuteman  Missile.  The  following  considerations 
were  appropriate  when  fabricating  the  experimental  pellets: 

1.  The  use  of  a solid  cylindrical  fuel  pellet  was  appropriate 
for  determination  of  data  necessary  to  achieve  the  first 
two  objectives  of  this  project.  The  use  of  the  actual 
Minuteman  fuel  pellet  geometry  would  be  necessary  to 
verify  the  final  mathematical  formulation. 
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2. 


Special  precautions  were  necessary  to  insure  a proper 
medium  was  available  to  collect  the  appropriate  data. 

The  actual  plasticizer  that  was  employed  in  the 
Minuteman  fuel  pellet  was  selected  for  inclusion  in 
the  experimental  pellet  while  the  other  ingredients  were 
carefully  selected  to  be  as  close  a facsimile  to  the 
real  fuel  pellet  as  possible(l). 

The  experimental  samples  were  constructed  with  the  following 
composition: 


Prepolymer  (R45M)  157^ 
Curing  Agent  (DDI)  37c. 
Plasticizer  (DOA)  67, 
Filler  (Ammonium  Sulfate)  76% 


These  constituents  were  stirred  until  a homogeneous 
mixture  was  attained.  Since  the  material  was  extremely 

viscous,  it  was  necessary  to  use  a vacuum  apparatus  to  remove  excess 
air  bubbles  in  the  mixture.  The  mixture  was  then  cured  at  60°  C 
for  24  hours. 

Obiective  Number  One  Consideration 

Experimental  data  collected  by  Jendrek  using  a chemical  mix 
similar  to  that  used  in  an  actual  fuel  pellet  were  available. 
Initially,  weight  loss  measurements  were  performed  to  determine  which 
mechanism  for  plasticizer  disappearance  seemed  appropriate  for 
further  study.  At  this  point,  three  principal  mechanisms  had  been 
theorized.  These  were  as  follows: 

1.  A diffusion  controlled  mechanism  wherein  only  the  physical 
process  of  plasticizer  diffusion  from  within  the  fuel 
pellet  to  the  surface  of  the  fuel  pellet  occurs. 

2.  A reaction  controlled  mechanism  wherein  the  plasticizer 
chemically  degraded  resulting  in  further  degradation  of 
the  fuel  pellet. 

3.  A combination  of  both  diffusion  and  chemical  degradation 
resulting  in  a net  loss  of  plasticizer  from  both  processes. 

Thus,  to  proceed  with  the  other  objectives  outlined  earlier,  it  became 
necessary  to  determine  which  mechanism  was  appropriate  for  further 


Mechanism  number  two,  a reaction  controlled  phenomena  was  first 
eliminated  for  several  reasons.  It  was  thought  that  the  possibility 
of  chemical  degradation  was  remote  since'  dloctyl  adipate  is  a 
relatively  stable  high  molecular  weight  organic  material.  The  low 
temperature  range  and  relatively  inactive  Ingredients  used  in  the 
fabrication  of  the  fuel  pellets  were  not  conducive  to  producing  any 
type  of  chemical  degradation.  In  addition,  weight  loss  measurements 
indicated  that  material  was  disappearing  from  the  experimental  fuel 
pellets.  Since  a reaction  controlled  mechanism  would  not  yield  a 
gross  weight  loss  unless  the  reaction  products  themselves  diffused 

' through  the  material,  it  was  again  implied  that  the  reaction  controlled 

mechanism  was  an  unlikely  candidate  for  further  study. 

* Mechanism  number  three,  a combination  of  diffusion  and  chemical 
degradation,  vas  also  eliminated  because  of  weight  loss  considerations. 
This  was  due  to  the  conclusion  that  a relatively  large  percentage  of 
dioctyl  adipate  was  disappearing.  This  could  be  explained  only  as  a 
diffusion  controlled  phenomena  and  not  a reaction  controlled  phenomena 
wherein  no  gross  weight  loss  would  be  encountered. 

Consequently,  it  was  determined  in  a subjective  fashion  to  pursue  the 
diffusion  controlled  mechanism  as  a means  to  describe  the  phenomena 
occurring  within  the  fuel  pellet. 

Obiectlve  Number  Two  Considerations 

In  order  to  eventually  produce  a useful  correlation  of  the  following 
form, 

Cs  = f (t,  r,  z)  (1) 

three  physical  parameters  describing  the  diffusional  process  appear 
in  any  attempt  to  mathematically  model  this  phenomena.  These 
parameters  are  Dj  (the  molecular  diffusivity  of  the  plasticizer  in 
the  solid  fuel  pellet),  Djj  (the  molecular  diffusivity  of  the 
plasticizer  in  the  gaseous  medium  surrounding  the  open  surfaces  of 
the  fuel  pellet),  and  M (a  Henry's  Law  type  distribution  coefficient 
correlating  the  concentration  of  plasticizer  across  the  surface  of 

* the  fuel  pellet).  As  very  little  information  was  available  on  the 
three  parameters  noted,  the  use  of  mathematical  modeling  techniques 
to  design  experimental  methods  for  determination  of  typical  values 
for  each  parameter  was  accomplished. 

The  physical  structure  dictated  by  the  experimental  design  and  used  to 
derive  the  solid  state  molecular  dlf fuslvities  was  a cylindrical 
pellet  surrounded  by  a plastic  sheath  with  one  end  coated  to  insure 
diffusion  in  only  one  axial  direction.  The  pellet  was  kept  in  an 
inert  nitrogen  atmosphere  and  periodically  weighed  for  gross  weight 
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loss  These  circumstances  were  forced  by  various  boundary 

condition  constraints  imposed  during  the  mathematical  modeling  of 
the  diffusional  phenomena. 


The  mathematical  model  and  the  techniques  used  to  develop  the  desired 
correlation  for  the  determination  of  Dj  are  contained  in  Appendix  A (2)  . 
The  correlation  itself  is  as  follows: 


< -2  (-1)^  e-  {!T(2i  - 1V2]2 

\ ^ ^ g2i-l)  /2^Tr 


Cos 


TT  rL 
(2) 


where 


I>I  t 

17" 


Cs 

Cso 


z 

L 


(3) 


(Consult  the  Nomenclature  Table  for  definitions  of 
i various  symbols) 

I 

\ Mathematical  manipulation  of  Equation  (2)  yielded  a similar  expression 

I with  the  average  concentration  of  the  plasticizer  replacing  the  point 

* concentration  (Cs) . 

I 

The  experimental  data  taken  from  the  fuel  pellets  thus  allowed  determl- 

I nation  of  all  physical  parameters  in  Equations  (2)  and  (3)  except  Dj. 

? 

) 

While  the  previous  design  of  experiment  techniques  centered  on 
I developing  data  and  correlations  that  permitted  the  calculation  of 

? typical  values  of  Dj,  another  set  of  experimental  data  and  mathematical 

correlations  were  derived  to  permit  the  calculation  of  typical  values 
: of  M and  Djj. 

The  simplest  experimental  apparatus  that  permits  fairly  exact 
mathematical  description  and  results  in  a correlation  that  permits 
calculation  of  the  distribution  coefficient  (M)  was  constructed  by 
placing  a fuel  pellet  of  length  L in  a glass  tube  of  the  same  diameter 
as  the  pellet.  One  end  of  the  pellet  was  sealed  to  restrict  diffusion. 
Again  the  other  end  was  left  open  to  permit  plasticizer  diffusion  only 
in  the  axial  direction.  The  glass  tube  was  constructed  to  be  long 
enough  to  assure  no  plasticizer  would  diffuse  the  entire  length  of  the 
tube  after  leaving  the  open  surface  of  the  fuel  pellet.  As  in  the 
previous  model,  these  physical  circumstances  were  forced  by  various 
boundary  condition  constraints  Imposed  during  the  mathematical 
modeling. 
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The  actual  mathematical  mv/iel  and  techniques  used  to  develop  che 
desired  correlation  for  determination  of  M are  contained  in 
Appendix  B (3)  . The  final  correlation  is: 


r 


Csavg  — Cso 
- CsoM 


M + 


(1  + 


1 


L.  Jif 


(4) 


By  taking  weight  loss  measurements  and  measuring  various  physical 
characteristics  of  the  fuel  pellets,  it  was  possible  to  determine 
all  terms  in  Equation  (4)  except  M and  Dij.  Dj^  was  calculated  from 
an  empirical  equation  developed  by  Hirschfelder , Bird  and  Spotz  (4)  _ 
Consequently,  M was  determined  as  the  only  unknown  in  Equation  (4). 


RESULTS 

To  demonstrate  the  usefulness  of  Equation  (2) , values  of  Dj  have  been 
derived  from  the  experimental  data  available.  Table  1 shows  a cross 
section  of  the  values  obtained.  These  values  generally  agree  with 
published  values  of  solid  state  molecular  dif fusivities  * Values 
of  an  order  of  magnitude  of  lO"^  m^/hr  to  10"^  m^/hr  are  commonly 
reported  for  solid  .••.ystems  similar  to  the  system  presently  under  study. 

TABLE  1 


Equation  (4). 


Sample 

Time  (hr) 

Dj  (m^/hr) 

1 

1660 

9.1  X 10‘® 

2 

1870 

9.7  X 10‘® 

3 

2037 

9.8  X 10“® 

4 

2306 

9.6  X 10‘® 

distribution  coefficient  (M) 

were  calculated  using 

These  values  and  those  of  D 

TT  calculated  from  the 

Bird  and 

Spotz  correlation  are  shown  in  Table  2 

TABLE  2 

Sample 

Dfi  (m^/hr) 

M X 10^ 

1 

.194 

2.1 

2 

.194 

3.4 

3 

.194 

2.5 

4 

.194 

3.0 

(6) 


As  previously  stated,  the  values  presented  in  Tables  1 and  2 are 
presented  only  to  demonstrate  the  usefulness  of  the  correlations 
developed  earlier. 
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In  order  to  partially  verify  the  conclusion  deduced  earlier  wherein 
a diffusion  controlled  mathematical  model  would  be  the  most  appropriate 
as  a description  of  the  actual  phenomena,  Equation  (2)  and  the  average 
value  of  Dx  were  then  used  to  predict  gross  weight  loss  from 
independent  fuel  pellet  samples.  Figures  }.  and  2 indicated  that  a 
diffusion  controlled  model  is  adequate  to  approximate  the  phenomena 
that  occurs  within  a fuel  pellet  . 


FUTURE  CONSIDERATIONS 


While  the  work  necessary  to  accomplish  the  first  two  objectives  of  this 
project  are  essentially  complete,  the  work  to  complete  the  third 
objective  is  in  the  preliminary  stages  only.  The  most  feasible  model 
to  describe  an  approximation  of  the  actual  fuel  pellet  geometry  is 
briefly  discussed  in  Appendix  C.  Further  considerations  of  this  model 
were  postponed  until  completion  of  the  first  two  objectives  in  order 
to  determine  whether  the  relative  magnitudes  of  Dx  and  Dxx  would  pernit 
elimination  of  several  terms  as  insignificant  when  these  terms  are 
compared  to  more  significant  terms.  The  following  work  is  proposed  to 
adequately  verify  the  usefulness  of  any  correlation  derived  from  the 
model  suggested  in  Appendix  C. 

1.  Complete  the  mathematical  solution  of  the  proposed  model 
in  lieu  of  various  simplifications  that  may  be  available 
from  consideration  of  the  relative  sizes  of  Dx  and  Dxx« 

2.  Construct  several  samples  to  approximate  the  geometry  of 
an  actual  fuel  pellet  and  take  the  appropriate  weight  loss 
data. 

3.  Construct  several  standard  solid  cylindrical  fuel  pellets 
of  the  exact  same  mix  as  was  used  to  construct  the  pellet 
described  in  2 above  to  be  used  with  Equations  (2)  and  (4) 
for  the  calculation  of  Dx  and  M. 

4.  Use  the  correlation  developed  in  1 above  and  the  values  of 
Dx  and  M calculated  in  3 to  compare  to  the  actual  weight 
loss  curve  found  from  3 thus  verifying  the  usefulness  of 
said  correlation  to  accurately  predict  plasticizer  loss 
within  an  actual  fuel  pellet. 

Objective  number  four  can  then  be  accomplished  by  subjectively 
evaluating  the  terms  contained  in  the  prediction  correlation. 
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FIGURE  1.  Plasticizer  C 
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APPENDIX  A 


Derivation  of  a Correlation  for  Calculation  of  Solid  State  Molecular 
Dlf fusivlties  in  a Fuel  Pellet. 


Solid  Pellet 


Figure  3 

Assuming  the  concentration  of  the  plasticizer  in  a solid  cylindrical 
fuel  pellet  is  not  a function  of  radial  direction,  but  will  diffuse 
uniformly  in  the  axial  direction,  Fick's  2nd  Law  of  Diffusion  yields; 


Rearranging  Equation  (lA)  using  the  following  dimensionless  variables 
yields: 

^ Y = ^ ^ Y 

(2A.) 


where 


Applying  separation  of  variables  techniques  to  Equation  (lA)  results  in 
the  following  expression; 


Y » Cl  (C2  Sino<,fV,+  C3  Cosolru) 


At  z = 0 (»V,=  0),  the  flux  of  material  is  physically  restricted  to  zero; 
therefore 


Since  the  molecular  dlffuslvlty  in  the  inert  nitrogen  atmosphere  was 
several  orders  of  magnitude  greater  than  the  molecular  diffusivity  in 
the  solid  phase,  material  can  be  carried  away  from  the  surface  faster 
than  it  can  be  supplied  from  the  solid  core.  Therefore,  at  z » L,  the 
assumption  was  made  that  Cs  equals  zero.  Consequently  the  following 
dimensionless  boundary  condition  was  applicable. 

at  'A.  - 1,  Y = 0 /fi.x 
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The  Initial  condition  for  the  time  variable  is  as  follows: 


I 


I - 0,  Y 


Use  of  Equations  (5A) , (6A) , and  (7A)to  evaluate  the  contents  in 
Equation  (4A)  yielded  the  desired  result. 


Y . ^ - 2 (-  1)^ 

1-1  [[(21  - i)/23Tr 


^(2i  - 


X Cos 


(“r*) 


TT  IV 


Solutions  of  this  nature  are  generally  stable  after  evaluating  the  first 
several  terms  of  the  expanded  equation.  Similar  solutions  appear  in 
engineering  literature  for  various  physical  geometries  . 


APPENDIX  B 


Derivation  of  a Correlation  for  Calculation  of  Molecular 
Distribution  Coefficients 


Sealed 

End 


Fuel 

Pellet 


Stagnant  Air  Atmosphere 


Solid-Gas  Interface 


Figure  4 


Since  the  fuel  pellet  was  fabricated  in  a cylindrical  glass  shell, 
plasticizer  diffusion  was  permitted  only  in  the  axial  direction. 
Therefore,  the  concentration  of  plasticizer  in  the  fuel  pellet  was 
not  a function  of  radial  distance.  Applying  Fick's  2nd  Law  of 
Diffusion  to  the  physical  system  described  above  results  in  the 
following  expression; 


(IB) 


The  following  initial  boundary  condition  applies  for  the  solid 
phase: 


at  t *■  0,  Cs  • Cso 

(2B) 

The 

axial  boundary  conditions  are: 

at  z = - , Cs  =>:  Cso 

(3B) 

at  z = 0 , Di  _ jj 

dz 

dCg 

dz 

(4B) 

Similar  equations  can  be  written  for  the  mass 
the  gaseous  side  of  the  solid-gas  interface. 
2nd  Law  of  Diffusion  yields: 

transfer  phenomena  on 
Again  employing  Flck' 

^ Cg  _ Dtt 

^ 2 Cg 

r 

^ t 

• (SB) 

The  initial  boundary  condition  for  plasticizer  diffusion  in  the 
stagnant  gaseous  region  is  as  follows; 

at  t « 0 , Cg  « 0 (BB) 

Axial  boundary  conditions  are; 

at  z - , Cg  - 0 (7B) 
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i 


at  2*0 


Di  d Cs 
d z 


Dji  dCg 
dz 


(SB) 


Cg  » MCs 


(9B) 


Applying  the  appropriate  Laplace  Transformation  techniques  to  this 
series  of  equations  results  In  the  following  expressions: 


Cs  - Cso  ^ 1 + erf  (z/\j4DTt) 

“CsoM  M + (Dj/Djj)^ 


(lOB) 


and 

Cso  (1  - erf  (2/\j4DTt) 

l/M  + (Djj/Di)^  (IIB) 

r 3^ 

These  expressions  are  similar  to  those  predicted  In  other  studies 
for  a similar  system. 


Integration  of  Equation  (llB)  and  application  of  the  appropriate 
mathematical  techniques  results  In  the  following  correlation  for  the 
average  concentration  of  plasticizer  at  any  time: 


Csave  - Cso 
- CsoM 


1 

M + (Dj/Djj)^ 


(1  + 1/ j4Dit  + 1/Lw%) 

(12B) 
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APPENDIX  C 


Proposed  Multi-Region  Model  for  Diffusion  of  a Plasticizer  from 
a Fuel  Pellet 


SECTION  1 


SECTION  2 


Figure  5 

Equations  of  State  for  Each  Section: 


For  Section  1 

a. 


+ i 


Cs 


Cs 


Dl 


^Cs 


For  Section  2 


^*7  “ii  ^ 


For  Section  3 


II 


^Cn 

Tt 


Initial  Boundary  Conditions 
at  t = 0 


Cs  • Cso 
Cg  ■ Cn  = 0 


At  z 


dCg 

dz 


- L 
0 . 


dCs 

dz 


At  z 


Dt  dCs 

■*"  J _ 


dz 


+ Dll  ^ 

dz 


Oil  dCr^ 
dz 


Cn  - M Cs 


Cn  » Cg 
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SECTION  3 


(1C) 


(2C) 


(3C) 


(4C) 

(5C) 


(6C) 


(70 

(80 

(9C) 


Separation  of  Variables  and  Use  of  Various  Simplification  Techniques 
will  yield: 


1 

1 

1 

© s ■ 

f(Cs)  - Ks 

Cos(Uz) 

(Jo(^r)  - Jx(^R2) 
Y1OVR2) 

Yo(Xr)) 

(IOC) 

©g  - 

f(Cg)  - Kg  e-'»^Ilt 

Cos(U' z) 

Jo(Xr) 

(IIC) 

©n  - 

f(Cn)  - f(e^  , Cos  z , 

Jo(r)) 

(120 

where 

U,  - f(M,  Di,  Djx)  (13C) 


I 
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NOMENCLATURE 


1. 

Cs 

- 

concentration  of  plasticizer  In  the  fuel 

pellet 

2. 

t 

- 

time 

3. 

r 

- 

radial  distance 

4. 

z 

- 

axial  distance 

5. 

Dl 

- 

molecular  dlffuslvlty  of 

the  plasticizer 

material  In  the 

fuel  pellet 

6. 

Dll 

- 

molecular  dlffuslvlty  of 
stagnant  gas  atmosphere 

Che  plasticizer 

material  In  the 

7. 

M 

- 

distribution  coefficient 

8. 

1 

- 

dummy  variable 

9. 

Y 

- 

Cs/Cso 

10. 

'V 

- 

(Di  t)/L 

11. 

VU 

- 

z/L 

12.  L - length  of  the  fuel  pellet 

13.  Cso  - Initial  concentration  of  the  plasticizer  In  the  fuel  pellet 

14.  Csavg  - average  concentration  of  the  plasticizer  In  the  fuel  pellet 

15.  Cl,  C2,  C3,  ot,  - separation  constants 

16.  Cg  - concentration  of  plasticizer  In  the  stagnant  gas  atmosphere 

17.  Cn  - concentration  of  plasticizer  In  Section  3 

18.  Os  - Cs/Cso 

19.  Og  - (Cg  + Cso  M)/Cso  M 

20.  ©n  - (Cn  + Cso  M)/Cso  M 

■'  21.  u,  u',.B»^>X»  Ks»  Kg  ■ separation  constants 
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